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I. Introduction 

Intuitive discussions of phenomena in quantum field theory frequently take 
the form of describing them as properties of the vacuum. The vacuum, it is 
then argued, is just like any other medium, so we can use the intuition de- 
veloped in other fields of physics to get an understanding of quantum field 
theory. Perhaps the first example of this was the interpretation of unoccu- 
pied negative energy states in Dirac's relativistic electron theory as positrons. 
Today the obvious analogy is the behavior of electrons in semiconductors, 
though historically it seems that the concept of holes developed independently 
in quantum field theory and solid state physics |l|. Another example is the 
Casimir effect, which is explained in terms of the zero-point fluctuations of 
the quantized electromagnetic field. A third example is spontaneous symme- 
try breaking of the vacuum, where the intuition is supplied by the behavior 
of (anti-)ferromagnets and (for local symmetries) superconductors. As a fi- 
nal example we have the charge-screening properties of the vacuum. Already 
in the 30s it was realised that the QED vacuum behaves just like an ordi- 
nary dielectric medium |J: Charges placed in the vacuum get screened by 
polarization effects. Around 1980 Nielsen, ||, and Hughes, ||, showed that 
asymptotic freedom, the antiscreening property of the QCD vacuum, could 
be interpreted as (color-) paramagnetism.^ In this paper we will show that 
the calculation this result rests upon, the change in vacuum energy due to an 
external (abelian color-) magnetic field, can also be used to derive the QCD 
trace anomaly. The point is that the energy is all we need in order to ap- 
ply some elementary arguments borrowed from the electrodynamics of media 
to get the stress tensor, and hence the anomaly. The vacuum is indeed just 
another (color-) electromagnetic medium. 

Our paper is organized as follows. We first (section 2) review the medium 
picture of screening (QED) and antiscreening (QCD). In the course of this 
discussion we quote the the result for the vacuum energy in QCD and introduce 
the renormalization group equation for charges. In section 3 we then use the 
expression for the vacuum energy and some "elementary physics" arguments 
to deduce the the stress tensor for QCD and calculate the anomaly. In these 
two sections we have tried to keep the discussion at an elementary level and 
to explain those results that we have just quoted from the literature. Finally, 
in section 4, we conclude with some comments as to how our derivation of 
the trace anomaly reflects the breaking of scale invariance. In appendix A 
we give a standard derivation of the trace anomaly and compare some details 
with our derivation. We also briefly discuss higher order effects and some 
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other technicalities. In appendix B we justify an argument based on Lorentz 
invariance used in the text. Both appendices requires more knowledge of 
quantum field theory than the main text. 

Although this work is a rederivation of a known result, we hope that it 
has a value beyond the purely pedagogical. The vacuum as a medium has 
been a fruitful picture in trying to develop intuition about non Abelian gauge 
theories. In order to make good use of the medium picture, and avoid pitfalls 
due to false analogies, we think it is of value to establish its limits as precisely 
as possible. 

2. Charge (anti) screening 

As a preliminary to the dicussion of the trace anomaly we wish to remind 
the reader about the intuitive interpretation of charge (anti-)screening. Start- 
ing with QED, we consider two heavy particles of charge e, separated by a 
distance r. Classically the potential energy of the pair is given by 



E(r) 
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where a = e 2 /4ir. This result is modified in QED. To 0(a 2 ) the potential 
energy is instead given by 
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where we have introduced an ultraviolet cut-off Q. The usual charge renor- 
malization in QED is obtained by taking 
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as the physical charge. This amounts to removing the cut-off dependent term 
from equation (3), and in the limit r — > oo, we get 
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The "effective charge" defined by 

7-1/ \ a eff 
E(r) = — — 
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depends on the distance. The intuitive explanation of this is that the QED 
vacuum behaves like a polarizable medium. In Dirac's original picture the 
vacuum has a homogeneous distribution of negative energy electrons. When 
an external charge is introduced, the charge distribution is modified. More 
negative charge accumulates around a positive external charge, less around a 
negative external charge. The effect is that the external charge is screened. 
This way of looking at the problem is not explicitly charge-symmetric, but the 
same intuitive picture works if we prefer to view the vacuum a consisting of 
(virtual) electron-positron pairs. 

The quantitative result for e(r) given in (|5]) has the property that e — > 1 as 
r — > oo. This limit is finite because the electron mass m/0 and we have used 
this to renormalize the electric charge in such a way that a means the charge 
at the "preferred" distance r = oo, but this is clearly a matter of choice. A 
situation more like QCD arises when m = 0. In the limit m <C 1/r, we get 
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Note that both terms are logarithmically divergent as m —>■ 0, but that the 
sum is finite. The charge renormalization can now be done in the following 
way: First, note that in the limit m — > 0, 

2a 2 r 
r8E(r) - r 5E(r ) = -— In - , (8) 

so we can write the total potential energy as, 

E = - + 8E(r) 

r 
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Then we define 



to get, to 0(a 2 ), 



a(r ) =a + r SE(r ) , (10) 



£(r) = ^(l-^Mln^) . (11) 

It is no longer possible to take the limit r — > oo in defining the charge, instead 
we pick an arbitrary scale r . The dependence of the charge parameter on 
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the defining distance is governed by the statement that the energy E(r) is a 
physical quantity and should be independent of arbitrary choices like r . Thus 

(12) 



dE 
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which, to this order in a, gives 
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This is the renormalization group equation for a. It is conventional to write 
renormalization group equations in terms of derivatives with respect to a mass 
scale fj, rather than a length scale. This merely reverses the sign of the right 
hand side of the equation and we get 

da n . . 2a 2 .„ ,. 

"Tr e(a)= ^ ■ (14) 

The description of screening, as formulated above, is in terms of electric 
fields from a point charge. One can also discuss the effect in terms of the 
response of the vacuum to an external homogeneous electromagnetic field. 
Since an external electric field always gives a non- vanishing probability for 
pair creation, it is more convenient to discuss the response to an external 
magnetic field and calculate the magnetic permeability fi (not to be confused 
with the mass scale in the renormalization group equation!) rather than e. We 
can then deduce e from \x if we assume that the vacuum is Lorentz- invariant. 
If that is the case (for further discussion, see Appendix B), then we should 
have 

€fi=l , (15) 

so, if /i > 1 (paramagnetic vacuum), then e < 1 (antiscreening) and if /x < 
1 (diamagnetic vacuum), then e > 1 (screening). The strategy is thus to 
calculate the energy of the vacuum in an external magnetic field and read off 
/i. This was done in ||, ||, for the case of a constant homogeneous color 
magnetic field belonging to an abelian subgroup of the color group, i.e. of the 
type 

A a {x) = A(x)T a (16) 

B a (x) = V x A(x)T a (17) 

where T a is a constant color matrix which is usually taken to be diagonal (i.e. 
for SU(2) proportional to a 3 and for SU(2) e.g. proportional to A 8 ). We refer 



the reader to the above references for the details, here we will just outline the 
calculation and give the final result. 

The vacuum energy is calculated by summing the zero-point energies for 
the modes of all fields. For the gluon field this is ^ for each mode, the usual 
harmonic oscillator zero-point energy, for the quarks the exclusion principle 
(technically, the fact that fermion creation/annihilation operators satisfy an- 
ticommutaton relations instead of commutation relations) reverses the sign of 
this, so the contribution is — ^ for each mode, counting particles and antipar- 
ticles separately. In the Dirac sea picture, one can also think of the fermionic 
part as coming from the filled negative energy states, where each mode con- 
tributes — Two. 

The frequencies for the various modes are essentially the energy eigenvalues 
for particles with spin 1 (gluons) or spin 1/2 (quarks) in an external magnetic 
field. As is well known from elementary solid state physics (in the case of 
spin 1/2), both the spin and the orbital magnetic moment give contributions. 
Remember that without spin, an ordinary electron gas would be diamagnetic, 
with the spin included it is actually paramagnetic. In the case of the vacuum 
medium the situation is reversed. The vacuum electrons give a diamagnetic 
contribution due to the extra sign coming from the exclusion principle. For 
gluons the story is similar, but since there is no extra minus sign, their contri- 
bution is paramagnetic. Of course, when we sum over all modes, the energy 
actually diverges. One has to introduce a cut-off fi and subtract a piece oc fi 4 
that is independent of the external field. The result for the field- dependent 
piece in QCD with Np species of mass-less colored fermions is 
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where g is the QCD coupling constant and V is the volume of the box whose 
energy we have calculated. The first term in ( Jl8| ) is the classical contribution 
to the energy, just like the Coulomb contribution in the case of external point 
charges. We still have an explicit fi in this expression. We can handle it in a 
way that is similar to what we did in the point charge case above and which 
allows us to let Q — > oo. By rearranging the expression in ( 15 ) we can write 
it as follows: 
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We now renormalize the field strength, in analogy with (10), by defining 
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where all field strengths are renormalized ones. Just as we did for the charge 
in QED, we can now derive a renormalization group equation for B(gB ). The 
result is 
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If we remember that gB is a physical quantity, we can readily convert this to 
a renormalization group equation for g: g r B r = gB, so 



gB, 
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where the coupling constants are renormalized ones. We have introduced the 
conventional /3-function, and the factor 1/2 in the last equality arises because 
gB Q has dimension (mass) 2 , and, as we said above, it is conventional to use 
the derivative with respect to something with dimension mass. 

To deduce a permeability from the vacuum energy, one usually defines // 
by 



u = -HB = —B 2 

2 2/i 



(24) 



where the first equation has been taken from the electrodynamics of linear 
media. This is strictly speaking not quite correct, and we will return to the 
relation between H, B and u in the next section, but for the present purposes 
it suffices. We get 
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so we see that if Np < 16, \x > 1, the QCD vacuum is paramagnetic, and, by 
the previous argument e < 1, the QCD vacuum antiscreens color charge. 



3. The trace anomaly 

The point of view we wish to emphasize in this paper is that the QCD 
vacuum is like any other nonlinear medium. In particular, once we have the 
energy as a function of the appropriate variables, we can use elementary vir- 
tural work arguments to deduce the stress tensor and thus calculate the trace 
anomaly. However, before turning to that argument we shall ask what the 
medium we are talking about really is. 

In ordinary electrodynamics for media, the inhomogeneous Maxwells equa- 
tions are given in terms of the fields D and H and the corresponding macro- 
scopic currents and charges. The induced current and charge distributions in 
the medium is taken into account via polarization and magnetization which 
relates the D and H fields to the fields E and B. Since the fields and the 
macroscopic currents are related by Maxwells equations, we can choose to ei- 
ther consider the currents or the fields as fixed. When discussing the vacuum 
the situation is less clear, since there is no obvious separation between external 
(or macroscopic) and internal (or microscopic) currents. Roughly speaking, it 
is the virtual particles that constitute the medium, but we have to make this 
statement a bit more precise. In QED the situation is rather simple if we con- 
sider the case of an external field and work only to lowest order in the coupling 
constant. In this case the only virtual particles are electron - positron pairs, 
and the calculation is exactly the one given in the previous section. Since 
there are no virtual photons, there is no ambiguity in how to define the back- 
ground field. Note that if we had choosen to start from external currents of 
electrons rather than from a background field, we would have had to face the 
problem how to distinguish the "real" (on shell) electrons in the currents from 
the virtual (off shell) electrons in the vacuum medium. 

In the case of QCD the situation is more subtle. Here there are virtual 
gluons produced to lowest order, and these have to be distinguished from the 
gluons that constitute the external field. This is done by employing the so 
called background field method. (For a review see e.g. 0). This allows for a 
separation of the gluon field in a background and a quantum part. Only the 
quantum part appears in loops and needs to be gauge fixed. The background 
(or external) field is treated as classical, but due to the quantum loops it 
gets a non-quadratic (but gauge- invariant) effective action. The separation 
between background and quantum part of the gluon field depends on the gauge- 
fixing and renormalization procedure and is thus to a certain extent arbitrary. 
Nevertheless it will provide a strict definition of macroscopic and microscopic 
and thus allows us to use results from the electrodynamics of media, just as 
in the case of QED. 

In order to identify the background field, which was introduced in (17), and 
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that appears in the expression (21) for the vacuum energy, we first note that the 
field (17) which enters our expressions is the curl of the vector potential (16). 
Thus it is divergence free, and should be a B rather than a H-field. Second, 
our complete hamiltonian density, including the background field energy, is 
(for the time being we put a tilde on the field to remind us that our aim is to 
identify it) 

B 2 

71 = -y + ^matter (P ~ 9 Equant ~ ffA) , (26) 

where 7i m atter is the hamiltonian density for quarks and (the quadratic part of 
the) hamiltonian for gluons in an external magnetic field. For a constant and 
homogeneous external field we can choose the vector potential A = |B x r. 
Varying the field we get 

Su vac = 5(vac\H\vac) = (B - M) • SB , (27) 

where we have identified the magnetization, 

Ml I ''-matter i \ / i *■ -i \ /oo\ 

= (vac\ rx = — \vac) = (vac \- r x i ivac) , (28) 

\ I 2 gA i / \ i 2 ji / v i 

where j is the current operator. In the non relativistic case the right hand 
side in (28) reduces to a sum of the well know expressions for the orbital and 
and spin magnetic moment densities. In the relativistic case, relevant for light 
quarks, there is no clean separation between orbital an spin contributions on 
the operator level, but the distinction can still be made by considering the 
following expression for the (unrenormalized) one loop vaccum energy 0] 



1 1 a 2 s / 1 \ U 2 
l _B 2 --B 2 (-l) 2a ^ Y m 2 --)l n| — | , (29) 
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which is valid for spin s — 0, |, 1, and where g depends on the color repre- 
sentation of the particles. The term ^ comes from the orbital motion and 
is common for both scalar, spinor and vector particles, while the term m 2 , 
where m is the component of the spin in the direction of the B field, does 
depend on the spin. For gluons this expression coincides with (18) if g 2 = \g 2 -, 
corresponding to adjoint gluons. 

To identify H and B we compare Su vac with the work done as the magnetic 
field is varied. As the field changes an electric field E is induced and the work 
done on the system is M, 



SW = -St f d 3 xE ■ j = f d 3 x 5 A ■ j 



= f d 3 x5A-(Vx!I)= f d 3 xH-5B , (30) 

where we have used Maxwells equation 

VxH = -j , (31) 

c 

and j is the current density that serves as external source for the magnetic 
field. The work per unit volume (30) should equal the change in energy (27): 

H 5B = (B - M) SB . (32) 

Since B = H + M, we identify 

B = B , (33) 

and 

The work we discussed above is the total work on both medium and sources of 
the magnetic field. To be able to deduce the stress tensor we need to consider 
a situation where we can divide this work into work done on the medium and 
on the sources. This is what we are going to do next. 

In order to discuss the work done in changing the magnetic field it is conve- 
nient to have a specific geometry in mind. We consider a cylindrical solenoid 
of length L, radius R and with N turns. A current I through the solenoid 
generates a H-field, which, in the limit L — > oo, N/L finite, is homogeneous, 
parallell to the axis of the solenoid and has magnitude 

liV 
H=--I . (35) 

c L 

This solenoid encloses our medium, i.e., the vacuum. We can define the stress 
tensor in terms of the virtual work done as we change the enclosed region while 
keeping the current constant. (The more familiar electrostatic analogue of this 
procedure is to consider a region bounded by capacitor plates and the work 
done as the plates are moved while keeping them at constant potential 0.) 
First, it is clear from the symmetry of the problem that the only non- vanishing 
components of the stress-tensor Ty are T zz (where the z-axis is the axis of the 
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solenoid) and T xx = T yy . Thus Ty simply acts like an anisotropic pressure. As 
we change the solenoid we have to do a virtual work per unit volume 

owi = — p— — =—l zz — i X x ^r~ ■ <jo) 

In addition to this work on the medium we also have to do some work to keep 
the current constant. As the solenoid is modified the magnetic flux changes 
and an emf, 

« = -!>*■) , (3?) 

is induced. This emf does some work which has to be compensated for to keep 
the current constant. The additional work per unit volume is given by ([§], p 

214) 

' 5(-nR 2 LHB) , (38) 



nR 2 L 

where we have used that / is constant and the relation between / and H. The 
total work per unit volume is thus 

6w = 5wi + 5w 2 = (HB-T zz ) 6 -^ + (HB-Txx)^ + 5(HB) 

, TT2 dB m JL /ttt ^ m S 28R 
= (-H 2 —-T ZZ ) T + (HB-T XX )— , (39) 

where we have used that H, and consequently B, is independent of R and also 
that H oc L' 1 , i.e. L5H = —H5L). This work should equal the change in 
energy of the system per unit volume. This change has two terms. The first is 
due to that the solenoid encloses a different volume, the second is due to the 
changed magnetic field. The result is 

5V s 25R . TT2 dB^6L 

5 w = u— + 5u = u^- + {u-H 2 —)— . (40) 



In the last step we used that 5u = H5B, i.e. , equation (p4|), and again that 
H oc L _1 . We can conclude that, 



T T: -u 

r 
- yy 



zz 

T xx = Tyy = HB-u , (41) 
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and hence the trace of the stress-energy tensor is 

T» „ = Too - T xx - T yy - T zz = Au- 2HB . (42) 

From section 2 we know that 

B 2 33-2N F g 2 B\ , gB . 

u= — + ^-- n— , 43 

2 48tt 2 2 l gB l ' v ; 

and thus the final result 

3 3 - 2N F 2 2 /3(g) 

^ — Q d = 

48vr 2 2# 

where we have used the lowest order result for the /^-function 



T% = - AO _ 2 9 2 B 2 = ^F^F, V , (44) 



_, N 33-2iV F 3 
/3(5) = ^g^-^ 3 • (45) 

We have thus recovered the usual result for the trace anomaly. Note that it 
was important in this derivation to have the correct relation between u, B 
and H. If we had applied linear-medium relations, the anomaly would have 
vanished. The entire effect comes from the non-linearity of the vacuum. 

Readers wishing for a comparison of the X^- we derived with a conventional 
field theory version are referred to Appendix A, where we also show that the 
coefficient in front of F^ V F^ V is indeed related to the /9-function. 

4. The trace anomaly and scale braking 

In the standard derivation of the trace anomaly, the relation to sponta- 
neously broken scale invariance is emphazised. At a classical level the fields 
scale according to their canonical dimension, 

B(ar) -»• \- 2 B(\x) (46) 

E(x) -► \- 2 E(\x) (47) 

under the rescaling x^ 1 — > Xx^ of the coordinates. Since the classical La- 
grangian for the electromagnetic or gluon field does not contain any dimen- 
sionful parameter, this scale transformation is a symmetry and there is a corre- 
sponding conserved current, the so called dilatation (or scaling) current given 
by 

JSu. = T> U ■ (48) 

Since the energy-momentum tensor is conserved, we get 

dA. = T* , (49) 

12 



so if the dilatation current is conserved, the trace of the energy-momenum 
tensor vanishes. (For a more careful discussion of this statement, see fp~0f| - ) 
The trace anomaly thus implies that the dilatation current is not conserved. 
The classical scale invariance has been broken by quantum effects. How did 
this happen? A look at our calculations shows that we at some stage had to 
introduce a scale to be able to define the theory. In the case of external charges 
in QED we needed the lengh scale r$ to define the charge, and in the case of 
an external homogeneous magnetostatic field we introduced the scale in the 
guise of a reference magnetic field gB$. In both cases we made the point that 
this reference scale (tq or gB$) was arbitrary and we used this fact to derive 
renormalization group equations. However, despite the arbitrariness of the 
scale, it had to be introduced, and its presence destroys the scale invariance. 
As we change reference scale and coupling constant, there is a combination of 
them that remains unchanged: A renormalization group invariant scale that 
can be used to parametrize the theory. By integrating the renormalization 
group equation 

"djl = Pi9) ' (50) 

we can trade the arbitrary scale /i for a (dimensionful) integration constant A 

A = fie J 9< A > ^ . (51) 

One might ask what is the advantage of having replaced the scale \i by the 
(seemingly) equally arbitrary scale A. This becomes clear if we reexpress the 
vacuum energy u vac in (21) in terms of a A# defined via (51) with /i 2 = gB : 



3 3 - 2JVp A , 

1927T 2 



uvac = TT^r^K (52) 



We see that A B is directly related to a physical quantity, while \x is not. In fact, 
equation (51) describes a a curve in the g — \i plane, and all points on this curve 
corresponds to the same physics (masses, scattering amplitudes etc. ). To get 
different physics we have to change to a different curve, i. e. change the value of 
A. It is, however, important to realize that the actual numerical value (in e.g. 
MeV) of A depends on how the renormalization is performed. That is why we 
used a subscript B on the A in (52). There are many different prescriptions for 
defining A, but they can all be related, even though the required calculations 
sometimes can be lenghty.Q 

2 For instance, to relate the A that occurs naturally in lattice calculations to the one 
used in perturbative QCD, one must perform a rather complicated perturbative lattice 
calculation. 
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Classically the theory is scale invariant and is characterized by a dimen- 
sionless parameter, g. Quantum mechanics breaks the scale invariance and 
transmutes g to a dimensionful parameter, A! In QCD with mass less quarks 
this is the only parameter that sets the scale for hadron masses, crossections, 
etc., and it is used to parametrize theoretical predictions for QCD processes. 

To summarize: The trace anomaly implies that scale invariance is broken 
by quantum corrections, and this occurs because, in one way or another, the 
renormalization procedure forces us to introduce a scale to define the theory. 
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Appendix A. 

In section 3 of this paper we used elementary arguments to deduce the 
stress-tensor from the energy density. In this appendix we first show that 
this way of deducing the stress-tensor coincides with the conventional field 
theory definition in terms of a variation of the effective action with respect to 
the metric. By combining the effective action method with renormalization 
group arguments we then show that the coefficient in the trace anomaly is 
indeed related to the /3-function. (The argument here is the inverse of the one 



presented in [[11]]. We thank R.J. Hughes for informing us about this reference.) 
The starting point for our field theory considerations is the renormalized 
effective action 



(a.i; 



The stress-energy tensor is then given by[] 



2 5T dC 

T a p - —==-T-za - 2^—3 - g a pC . (A.2J 



-g 5g a P dg a P 

We are interested in T a p for a constant (abelian color) electromagnetic field, 
in which case C = C(x,y) where x = F a ®F a p and y = e al3lS F a/3 F^s / \/—g are 



3 Note that the introduction of a metric tensor is just a trick to calculate the symmetric 
energy-momentum tensor. With some care, the same result can be obtained by canonical 
methods. 
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the two possible invariant combinations of field strengths. We get 

dC dC 

T a p = 4F Q 7 ^>7"o- + 9a$ V-K 9a/3^ • (A.3) 

In the following we will ignore the dependence on y which is irrelevant for the 
discussion in this appendix. In a purely magnetic field 



and thus 





Foi 
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F- 

1 V 


= —€ijkB 


Too - 


= -C 




Toi = 
T - 


= 


~ B dB* ) + dB t Bi 



(A.4) 



(A.5) 

which we can compare with our "elementary physics" result from section 3 
(slightly generalized), 

r 00 = u 

(fll fill 

T tj = 5 l3 {-u + B k —)-—B r (A.6) 

Thus the results do agree if u = —C, as indeed it is [[|. We can now proceed 
to use the field theory calculation to show that the coefficient in the trace 
anomaly really is the /3-function. To this end it is convenient to rescale the 
gauge fields by gA^ — ► A^. Then the coupling constant only occurs as a 
coefficient in front of the gauge theory action, 

S = Jd A x{-±- 2 F^F, u + ...} , (A.7) 

and we can define an effective coupling constant in the renormalized effective 
action by writing 



/ 



d 4 xC =- I d^x —F^Fw . (A.8) 

A 9 2 e ff 



This expression defines an effective (or running) coupling constant g e ff which 
is a function of the original renormalized coupling g(fi), fi and the invariants x 
and y. Recall (section 2) that /x is the renormalization scale. Since according 
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to (A. 5) — C is nothing but the renormalized energy density of the vacuum, it 
is an observable and thus cannot depend on the renormalization point [i. 
Combining our equations we get 



rpa 



4 X g 9eff 

. 9eff dx 



-5-x . (A.9) 

A 9eff 



Dimensional analysis implies that the /i-dependence of g e ff is of the form 

X V 
9eff = 9eff(-4, ~i,9(lJ>)) , (A.10) 

and since C and thus g e ff does not depend on the arbitrary scale //, we get 
^ du ~ dx ^du da ' [ ' 



and thus 



J dfi dg 2g efJ 



r °'a = ^ a 03 F ^ ■ ( A - 12 ) 



The /3-function is actually the component of a vector field in the space of 
coupling constants: 



It follows that 



"!^ =«*»> • « a - i4 » 



and thus 



(Hgeff 

2 9eff 



T a a = ^f±F^F, u , (A. 15) 



which, up to the rescaling of the fields and to lowest non-trivial order in pertur- 
bation theory, agrees with the result derived in section 3 and moreover makes 
the relation to the /5-function manifest. 

The vacuum energy derivation of the trace anomaly given in section 3 makes 
use of an expression for the energy that is essentially a one-loop result, whereas 
the arguments in this appendix seems to involve no such approximation. In 
principle we could calculate the background field effective action to any desired 
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loop order, and it can still always be written in the form (A. 8). For this 
conclusion we have to keep in mind that there are two ways of presenting 
the anomaly. One is to establish an equality between the operator Tjf and 
the renormalized field operators of the theory. The other is to deal directly 
with matrix elements for T^ for states where there is an external gauge field. 
Both the vacuum-energy derivation and the background field effective action 
derivation follow the second way. It is known that as an operator equality 
(A. 15) is not true in general (beyond one loop). For further discussion, see the 
review fl2j by M. A. Shifman. 



Appendix B. 

In section 2 we appealed to Lorentz invariance to get a relation between e 
and \x. Actually it is not quite obvious that we can do so. After all, we want 
the relation in the presence of an external field so the physical situation is not 
Lorentz- invariant. We can use the effective action formalism to see under what 
circumstances we can still use the naive argument. It suffices for our purposes 
to consider the case of a slowly varying external abelian color electromagnetic 
field. Then again C = C(x, y 2 ) with x = 2(B 2 — E 2 ) and y = 4E • B. (Since y 



is a pseudoscalar it can't occur in odd powers.) We get (13| 



<9E dx dy 2 



dc ac-. . ac 

tt- = -4t^-B - 8y— 
<9B ox dy 

We see that if E • B = and if we define 



H = -^--4— B-%— E . (B.i; 



D = eE, B = /iH , (B.2) 

we do get 

e/i = 1 , (B.3) 

(assuming the derivatives do not become singular as y — > 0). In particular, we 
can use this relation for a homogeneous magnetostatic field. 
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